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Abstract 


The first and second-order density matrices for a uniform 
electron gas are discussed in the limiting cases of weedc coupling 
(usual gas parameter r_ -» O) and strong coupling (r_ ->00), 
Attention is focussed on the Important pair distribution function 
and on the momentum distribution. If the high density form of 
the momentum distribution as given by Daniel and Vosko (I96O) 
is adopted, then conclusions can be drawn regarding the meaning 
of a Fermi surface in a system of interacting particles. 

However, the applicability of perturbation theory to the 
calculation of the momentum distribution is thrown into some 
doubt by a calculation we have carried out on a soluble 
problem. Here we find that the perturbative answer is not 
correct, but unfortunately the problem, that of non-interacting 
electrons in a magnetic field, is very different physically 
from the Coulomb correlation case. Further work therefore 
remains to be done on this point. 

Variational forms of second-order density matrices are 
discussed, suid the Euler equations are obtained for one 
possible scheme based on localized orbitals. 

Finally, some progress on the non-uniform gas is briefly 
reported. 
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1. Introduction 

The problem of dealing with the Coulomh Interad-ions 
between conduction electrons In metala has proved yer^ formidable 
and no complete solution has, as yet, been fo\md. In recent 
years, beginning with the Important Investigation of Macke 
( 1950 ), attention has been focussed primarily on the use of 
perturbation theory. Unfortunately, the perturbation expansion 
does not prove useful In practice In the range of real metallic 
densities, which are characterized by the usual gas parameter 
r_ lying In the range 2 < r_ < 5*5* 

O B 

In this report, we consider chiefly methods based on a 
variational approach to the correlation problem, although 
considerable use Is made of results obtained In the limiting 
cases of weak (r. 0) and strong (r. eo) coupling. 

Certain very Important questions connected with the meaning 
of a Fermi surface In a system of strongly Interacting particles 
are not finally answered, but some light Is thrown on the 
general problem In sections 2 and 3 of the Report. 

While most of the work has been concerned with the uniform 
electron gas, some progress on the non-tmlfom problem le 
described In section 6. Here, a good deal of work remains to be 
carried out. 

We should mention that our earlier work on spin density 
waves has not been Included In this Report, because, after our 
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yarlatlenal ealoulatlons were eoapleted acre powerful aethods 
ware developed by Overhauaer and ear awn work is therefore» at 
most, of hlatorloal Intoreat. 
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2. Coulomb eorrtlatlona In a uniform elactron gaa 

We discuss first the role of electron Interactions In a 
uniform system. Throughout, ne shall focus attention on density 
matrices, and in sections 2.2 and 2.3 «e present results In the 
Sommerfeld model which are valid for the limiting cases of low 
and high density respectively. Some Inferences concerning the 
Perml surface are drawn In section 2.4, but no completely final 
result on the change In the momenttm distribution as the 
Interaction strength Is varied can be obtained, as some doubts 
exist as to the validity of perturbation theory in calculating 
this quantity. Some examination of this point is carried out 
in section 3 therefore. 
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2.1 The Hamiltonian and definitions of density aatrlcea for uniform 

gas 

We require basically to oaleolate the ground-atate ware 
fiinotion t(ziX 2 ....Xj|) for a ayatem of N electrona. where 
denotnn the space eo-ordlnatea r^ and spin co-ordinates 
from the SchrDdlnger equation 


H# 


et. 


( 2 . 1 . 1 ) 


We adopt aa the Hamiltonian of the Sommerfeld model: 

H = - i ^ ^ - 5^ YJ 

i i<d 


^ II 


dTo dTi 
(ro-ri| » 


( 2 . 1 . 2 ) 


where the last two terms In (2.1.2) refer to the electron-positive 
background Interaction and the background self-energy, 
respectively. We now define first and second-order splnless 
density matrices Y(ri' ** 1 ) «nd r(r,' ra' ri ra) by the equations 


Y(ri'ri) - ^ j ♦*(ri'OiXa...Xjj)t(r,o,Xa...Xjj)dOidXa...dXjj, (2.1.3) 


r(ri'ra'rira) » j ♦*(ri *0,ra'0aX3...Xjj;t(ri0ira0aX3...Xjj) 

dOi dOa dX3...dXj^. (2.1.4) 

We note now that the diagonal elements of y *nd r have direct 
physical Interpretations. Thus Y(ri ri) gives the particle 
density, which In the Sommerfeld model Is simply the constant 
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3/4iVg*» while r^PiPaPira) is effeetlrely the prohahllity of 
electron eeperation or the pair distribution function* Then* 
as liayer ( 1955 ) was the first to show* the energy per particle 
t/B Is glyen by 


where F(2^) » F(|r*-r|) Is the pair function* normalized to unity 
for large |r*-rl . 

Since we shall utilize the results later* we note that the 
usual Hartree-Fock solution based on a determinant of plane 
waves yields the forms: 

k-» 3i(k^lr*-r|) 

T(r'r) = -p- p - p cos p)/p» 


( 2 . 1 . 6 .) 


F(r*r) ■ 1 


g (- ji(k,lr'-r|)-)t 


whore h^ !• the magnltudo of the wave vector at the Feral eurfaeo 
and is related to the mean interparticle spacing by 

kfPg * * (2.1.7) 

Using atomic units* It follows from (2.1.5) end (2.1.6) that the 
energy per particle Is 

1-1% 
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The low-density form of the pair function has been dlaeuased 
earlier by March and Young (1959) using the electron lattice 
model of Wigner (1938)* We simply summarize the results by 
saying that as r. becomes very large, the electrons relative 
to a given particle we have singled out, and placed at the 
origin of co-ordinates sit on the sites of a body-centred cubic 
lattice. As r is then reduced somewhat, the electrons can be 
thought of as represented by harmonic oscillator functlona 



The results obtained by Maroh and Young (1959) for the pair 
function are represented in curves 1 and 2 of Fig. 1, for cases 
r^ s 100 and r^ s 4. For comparison, the Fermi hole result 
given in (2.1.6) is also shown in curve 3. We emphasize that 
while cxurve 1 should be reliable, curve 2 represents eui 
extrapolation beyend the range of validity of the low density 
form (2.2.1), and is given solely to show the qualitative Influence 
of varying r^. 

We shall now Indicate how this work may be generalized to 
yield the first-order density matrix and momentum distribution 
for a low-density gas. From the orbitals (2.2.1), centred on 
each lattice site, we build a Dirac density matrix. This will 
not of course, contain the condition of translational Invarlemce, 
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Fifi. I. Pair fuiictkim for an rleriron |ia» 

< urvr I. Low drii»ity form for r^ JOU 

<iur>r 2. Extrapolated low denAily form for r, 9 4 

Cur^e J. Hartref'Koi-k pair function of eqii. (2.^). rorrert in limit r, ■* 0 
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that la the diagonal eleaent will not he a eonatant. To obtain 
a result consistent with the Sonaerfeld nodel we auet awerage 
over all positions, and then the desired flrst'-order density 
matrix Is easily shown to he 

rir'r) « ^ exp J |r*-r|*y (2.2.2) 

The form (2.2.2) as Is easily seen, satisfies all the essential 
conditions, hut Its range of validity la, of course, restricted, 
as we shall discuss helow. Hevertheless, In conjunction with 
the high density results of $2.3, It may he used to draw soae 
interesting conclusions about the Feral surface (see $2.4). 

At this stage. It la enlightening to examine the aoaentua 
distribution corresponding to the first-order matrix (2.2.2), and 
this may he found as follows. We require the occupation numbers 
P(k) of plane wave states where V Is the volume of the 

metal and thus we write 

T(p'r) . ^ ^ P(k)e"^^’*‘' e^^***, (2.2.3) 

k 

or remembering that the density of states In k Is (l/8x^) V and 
using Bauer's expansion for a plane wave as a series of spherical 
waves: 

00 

Y(r'r) “ uJt I P(k) 4*k»dk. ^2.2.4) 

It Is convenient at this point to measure k In units of the Fermi 
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momentum tliat is we write K s k/k^, and then (2.2.4) heeemee 

k * ein k. K|r*-r| 

Y(r*r) = 5 ^ P(K) ^^^,7.;, K* dK. (2.2.5) 

Inverting this relation and using (2.2.2) for Y(r*r) we find 


P(K) * ^ exp K*j (2.2.6) 

s 

where we have eliminated k^ using (2.1.7). As we have remarked 
earlier, the range of validity of (2.2.6) is restricted, heeause 
the orbitals t en different lattice sites have been assumed 
orthogonal, whereas this is only rigorously true in the limit 

r -» 00 . However, a rough estimate of the range of validity 

8 

may be obtained by noting that the occupation numbers P(K) 
must always lie between 0 and 1. Since P(K) as given by (2.2«6) 
has its maximum at K « 0, we must have that 



(2.2.7) 


The kinetic energy per particle, in the approximation implied 
by (2.2.6), is easily shown to be 



in agreement with Wlgner (1938). 


( 2 . 2 . 8 ) 

Later treatments based on more 
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careful study of the dynaalca of the Tlhratlng electron lattice 
(Coldwell-Horsfall and Marad.udln, I960; Carr, 1961) change (2.2.8) 
somewhat, hut will not affeet the owerall validity of the 
argument presented ahove. We note finally for the low—density 
limit, that the momentum distribution la a classical 'Maxwell- 
Boltzmann* form: this Is In sharp contrast to the hlc^ density 
results to which we now turn. 

2.3 Pair function and first-order density matrix In hUdi-denelty 
limit 

Using the perturbation theory of Qell-Mann and Brueekner 
(1957) In suitable form, we shall new discuss the way In which 
the pair function and the first-order density matrix develop 
from the Hartree-Foek forme given In §2.1 ogioatlon (2.1.6). In 
order to do so. It Is very convenient to Introduce a generalization 
of dexislty matrices to Include also the time: tMs will greatly 
facilitate the setting up of useful equations from which the 
matrices can be determined. As a simple example, let us 
consider the Dirac first order density matrix for a system 
described by a single Slater detenilaant formed from one-body 
orbitals t^(r). This may be written 

Y(r'r) « ^ ♦l*(r')+j(if)» (2.3.1) 

1 

where the summation Is over the occupied energy states. However, 
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If we now introduce a generalized first-order density aatrlx 


*18 t ^ I B t 

Y(r't'rt) = ^ ♦^*(r')ti(r)e ^ e (2.3.2) 

1 

where are the occupied energy levels, then this is easily 
shown to satisfy the equation 

Hy = 1 , t ^ f (2.3.3) 

and this may he used to calculate y. But in fact, the above 
time-dependent matrix, is a special case of the one-particle 
Green's function used hy the field theorists in their approach 
to the many-body problem (see. for example. Klein and Prange. 
(1958)). 

Indeed, the convenient way to make the connection is now 
to start out from the definition of the Green's functions and 
the equations which they satisfy. The Green's functions or 
propagators are suitably defined matrix elements of Heisenberg 
field operators taken between exact eigenstates of the system. 
Following Klein and Prange we may write the one-particle and 
two-particle Green's functions as 

G(xx') - 1 < H tU(x) ♦*(x')| H > (2.3.4) 

and 

G(xiXaX 3 X*) » i* < H T| t(xi)t(xa)t*(x*)t*(x3)j H >, (2.3.5) 

where T is the time ordering operator of Wick. Then if we denote 
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the particle Interaction hy TCria) we have the fellowing equation 
for Q(x, x') 

7* + i 0(x,x') " ^ j d*r*T(r-r*)0(r*'trt,r*t'r't')“**(*“*') 

(2.3.6) 

and this reducea to (2.3.3) when the interaction is neglected. 
Also we hare 

0(xiXaX3X4) « 0(xiX3)0(xaX4) - 0(xiX4)0(xaX3) + 

+ J d*7i d*ya d*Zi d*Za 0(xiyj0(xaya)l(yiya*i*a)0(2i2aX3X4>. 

(2.3.7) 

Here I is a complicated interaction operator, which depends on 
the Green's functions and therefore to obtain exact solutions 
seems out of the question. Howewer, as we discuss below, in 
the high density limit it appears possible to make progress in 
the electron gas problem by inserting the first approximation 
for G(xiXaX 3 X 4 ) inside the integral, provided we approximate 
sufficiently carefully to the interaction operator. By now 
writing down the exoectation value of the Hamiltonian, we make 
the Indentlflcatlon between the Green's functions defined here 
and the density matrices of $2.1. and the results may be stated 
as follows: 

Lt G(r t r' t') = - ly(r' r) (2.3.8) 

t'-»t+ 

and 
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Lt 0(r t p' t Pi' t* Pa' t') “ “ 2r(ri' Pa' p')* (2*3«9) 

t'-»t+ 

Proceeding to the hlgh'-denslty gas ppohlem we make use of the 
work of Kanazawa and Watahe (1960), which shows that In the high 
density limit we must approximate I by 

I(yiya*<*a) » 6*(yi-ai) 6‘*(y8-»a) IV(yi-ya), (2.3.10) 

where the modified Interaction V Is given by 

V(yi-ya) ■ ▼(yi-ya) + j d^x d^y v(y,-x) Oo(x-y) 

Oo(y-x) V(y-ya). (2.3.11) 

Here Oq Is the first-order Green's function In the Hartree-Pock 
limit and the form of the effective Interaction has been 
discussed by Hubbard (1957). Then, In the Oell-Hann-Brueekner 
approximation, 

0(rtr'trt*r't'') » Oo(rtrt‘')Go(**'tr't") - i 0o(rtr't'')0o(r'trt") 

~ I dVidVa ®o(rtyi )a»(r'tya)V(y,-y8)Oo(yirt'')0*(ytr't*) 

^2.3.12) 

*ij dVi^^ya Oo(rtyi)Oo(r'tya)T(yi-yi)Oe(yir't‘')Go(y»rt''). 

We attempt at this stage to evalnate the pair function In 
the approximation In which we consider the first three terms on 
the rlfi^t-hand side of (2.3.12). The first two terms lead to 
the usual Hartree-Fook result given in $2.1, eg^atlQn (2.1.6), 
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and wa write oat only the third tern explicitly, called 
in what followB. It la aoat conTenlent to proceed by Fourier 
analyela, and to write 

dg^Crt r' t rt* r' t") - j *4 ^C^Cgt-t"). 

(2.3.13) 

If we Fourier tranafom V and Oo, th«i the reault aay be written 


9^(0.) ■ J ^ dk* de* dw Qo(kc)Oo(k- 4 ,e^) 


V(Qw)®o(k'e' )ao(k~ 4 e*-Hr). 
Introducing the function 


(2.3.14) 


00 

4o(qF) ■ i J d*k J de ao(4’«-ka4^w)Oo(ke) (2.3.15) 


aa In Oubola (1959) we aay simplify (2.3.14) to read 
, . 1 r Qo*(4.w)/**k-q* 


( 2 . 3 . 16 ) 


Numerical ewaluation of this correction to the Hartree-Fock pair 
ftaotton Is now In progress. 

The method just described may alao be applied to the first 
order Oreen's function to yield the firat-order density matrix. 
Results essentially equlTmlont to the moasntam distribution of 
Daniel and Voalte are then obtained* 
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2.U The ■•■nlng of a Feral siarfaea 

The feet that we ean obtain the aonentua dlatrlhutlon In 
two limiting eaeea has interesting implications. For, as the 
work of Daniel and Vosko shows, the discontinuity in P(K) at 
the Fermi surface for r. a 0 is not remowed hut only reduced 
as the density is lowered. While this result seems not to 
haye been established rigorously outside perturbation theory, 
there seems strong emplrloal STidence from the sharpness of 
the Fermi surface in real metals that such a discontinuity does 
in fact persist. However, our low density treatment r^ oo 
could equally be made the basis of a perturbation treatment, 
and no slon of the discontinuity would then occur. Thus, the 
evidence seems clear that as we follow the momentum distribution 
as it develops from the two limiting cases of small and large 
r. there must come a orltical coupling strength, or a critical 
density, at which the discontinuity in the momentum distribution 
is reduced to zero. For lower densities, it then appears that 
the concept of a Feral surface will no longer be useful. We 
believe at present that the critical value, r^ say, will lie 
outside the range of real metallic densities, or in other words 
that r^ will exceed 5*5* No quantitative evaluation has so far 
proved possible however. Questions also remain as to the nature 
and order of the 'transition' occurring at r^. We do not Azpeet 
that the pair function will undergo any marked changes at r^ and 
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prolia'bly, to obtain a clear plet\ire of the tranaitlon. It will 
eventually prove necessary to exaaine the third and higher- 
order density matrices. 

Finally, to give an indication of the change in the 
momentum distribution as the density is varied, we have plotted 
in curve 1 of Figure 2 results for P(K) as given by Daniel and 
Voslco for the high density values (r. | 2) and by (2.2.6) for 
the low density region. Curve 2 for r. 3 16 should be viewed 

B 

with some caution, but the result for r. s 100 (curve 3) should 

represent a good approximation. The great difficulty, as we 

have stressed, is to obtain a sufficiently good approximation 

to the ground state in the region 2 < r„ < 5*5 and to locate 

the critical value of r_ with precision. 
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3. Test of perturbation method for calculating electron monentim 
distribution (with D. Hilton) 

Sines the ahoYs work was conpleted, we have been eoneemed 
about the possibility that the Daniel and Vosko theory of the 
Bonentun distribution of a hlfi^ density electron gas may not be 
physically correct. Thus, as they In fact point out. the 
discontinuity In their BoaentuB distribution Bay arise froB 
perturbation theory, because It Is already present In the zeroth- 
order problen. 

In view of the Inportance of this question, we have sought 
a soluable problen In which perturbation methods can be avoided, 
and yet the electron BonentuB distribution can still be calculated. 
It seems to us that there Is a problem; not. unfortunately, closely 
related to the correlation problea; which we can solve exactly. 

In which a perturbation Is applied to an electron gas. This Is 
the case of non-interacting electrons In a magnetic field, and 
In view of the interest In the BomentuB distribution calculation 
by a non-perturbatlve method, we discuss It In some detail below. 

3»1 Bloch matrix and Bonentum distribution for non-intip p^et^nff 
electrons In r magnetic field 

As remarked previously. It Is fairly clear that a momentum 
distribution calculated by perturbation theory will contain a 
discontinuity If the unperturbed solution does so. It Is not. 
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howevery clear that the exact aolutlon will be dlecontlnuoua. 

In order to investigate the exact solution for the problem of 
free non-interacting electrons in the presence of a magnetic 
field, we note that the exact Bloch density matrix has been given 
by Sondhelmer and Wilson (1951) for this case for an infinite 
number of particles. It may be written in suitable units (the 
magnetic field being in the Z-dlrectlon) as 

C(r* £ B) = fO) exp[- ^ (x*y - j*x) - {(x-x*)® + (y-y')*i 

+ (a-z')*/UB] (3.1.1) 

where 


g(P) « H eoth HB and f(B) * — ^ r eosech HB. 

(*B)* 

It is easily seen that when the magnetic field H tends to zero, 
we recover the field free solution 

. -(r-rO*/4B 

* 


(*B) 


(Note that here the energy is 
expressed in Rydbergs). (3.1.2) 
We shall work with the quantity 


!)(» f) • j 0(J5 S' y) 4S' 




where 


ik'r -Ik.r* 

mm mm 


0(4 S' D e " " C(r £* B) dr dr' (3.1.4) 

Then by use of the relation pointed out by March and Murray (i960) 


\ 

\ I 
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1 r er E P) , , 

p(e' E O ^ j ^ (3.1.5; 

O—loo 

we see that 

D(k fi 

<i ()5 «) - I P(S S' - 2ir j “f— *e* (3.1.6) 

0>l«o 

We must remark that It la only in apeeial eaaea aa, for example, 
whAn the eleetrona are free, that Q (]( Z) la the momentum 
diayrlbution; we will return to thla point later. 

If we define the mixed matrix 
f ik.r' 

C(S E P) = J e " C(e'E P) *£' (3.1.7) 

it la eaally aeen that 


D(S P) = C(k r p) 

We perform the tranaform (3.1.7) by use of the formula 


(3.1.8) 


j . J I H, t > 0 ( 3 . 1 . 9 ) 

and find that 
o(s r 0) . pi 

( 3 . 1 . 10 ) 


or 


1 -PV -(k/+k,*)/g 

“ c o l h " ^ ® e * y . 


(3.1.11) 
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Now 


1 

i 


^ -4 p as H -♦ 0, and 


lia cosh H p a 1, so that 
H -»0 

lim D(k p) a 
H-40 


which is the correct answer. For we have for the momentum 

distribution in this case 

O+ioo 

* sir / e^^ dp. (3.1.12) 

O-ioo 


When > 2^ we close the contour by a semi-circle at 
infinity in the right half-plane. The integrand 
vanishes in this semi-circle, and since no 
singularities are enclosed the R(k K) vanishes. 

On the other hand if k* < Zt 1^1^* integrand vanishes 
in a semi-circle at infinity in the left-half plane. However 
the contour now encloses the singularity residue 1 at the origin. 
Thus 


«) - db: / 


O+ioo 


e“P^* ar 1 


O—loo 


k* < 5 
k“ > j; 


(3.1.13) 


the correct result for non-interacting electrons in the absence 
of a magnetic field. 

When the magnetic field is not zero, we hare 


t 

i I 
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0+1- -(kjj^+k *)/g0) 

-« 

O—loB 

o+i« 3(U-kj*-H) -(k^a+ky*)/gO) 

■5^/ - p(, ! 

O—iee 

Aftikin we close the contour hy a seal-circle at Infinity In the 
right half plane. The Integrand ranlahes on this seml-olrcle 
If 

- k a - H < 0. 
z 

Thus 

Q(k «) - 0 If kj* + H > j:. (3.1.15) 

This result Is at first sight a little surprising: we would 
expect the discontinuity to persist In the z-direotion» we 
might also expect it to remain In the same place, i.e. k^‘ » K 
when k^^ + ky‘ s 0, for the longitudinal motions of the electrons 
eu?e imaffected hy a magnetic field. However, the result may he 
explained hy supposing electrons originally moving wholly in 
the z-dlrection (that of the magnetic field) to Jump into states 
energetically more favourable emd for which there are non-zero 
components or momentum in directions transverse to the magnetic 
field. 

Evaluation of the Integral (3.1.14) does not seem possible 
when the equality in equation ( 3 . 1 . 15 ) does not hold, for there 
are essential singularities of the integrand on the imaginary axis. 

?.2 





We are thus compelled to proceed as follows. We put 

g » 51 - ic^a , H; T, - (k^* + y )/»• 


We now have 

1 HP ag(o+it) ^-Ti tan(t+lo) 

^ - A/ . T grov iu 

O—loo —00 

(3.1.16) 


and we will attempt to prove this function continuous as a function 
of Ti hy showing that its derivative with respect to t) exists and 
is finite everywhere. To do so we take the derivative under the 
Integral sign. It may be noted that we have already Interchanged 
orders of integration: we have taken a Foiurier transform before 
the Inverse Laplace transform. However, in that case, there is 
less apparent danger than here, where we seek to prove the 
derivative always bounded; for we must guarantee that if the 
transform of the derivative is bounded, then so is the derivative 
of the transform. We give an argument to justify our procedure 
which can also be used, with the necessary cheuges. to justify 
changes of order of Integration. 

We note that the relation 
00 

£ p (5 Ti 1 C(j5 11 ^) (3.1.17) 

can be made rigorous when C is a Dirlchlet series; when there 
is a finite number of particles, and the energy spectrum discrete. 
We have the functions of the form 
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p(5 T, Z) 


( 3 . 1 . 18 ) 


C(g Ti 8) 

H(li Z) 



(3.1.19) 


( 3 . 1 . 20 ) 


(The 8^ are the one-particle energy leYela). 

Now when there is a finite namher of particles, the allowed 
values of momenta (and so 2^ and t)) form a discrete set. We 
must take finite differences: 


p(€ Ti + h 5) - P(C T) + 5) 

Ap - -5- . (3.1.21) 


and ohvlously 


00 

I Ap e"^^ az = AC/p, 

Jo 


( 3 . 1 . 22 ) 


Now having passed to the limit of an infinite number of particles, 
we have a guasl-contlnuous function in the limit of the sets of 
functions obtained by Joining neighbouring allowed values of n by 
straight lines. For such a function, say P, the derivative is 
well approximated to be the finite difference of p for a finite but 
very large number of particles, l.e. if 

p^ s iim pjj, N the number of particles, 

N-ho 


XT = lia Ap„. 
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Such considerations must hold for both sides of equation (3*1 *22). 
Thus 




Our prciblem now consists of proving the boundedness of 


(3.1.23) 


We will use the second mean value theorem, which states that If 

^(x) Is a monotonlc function, and f(x) any other function, then 
b £ b 

j ^(x)f(x)dx = ^(a) j f(x)dx + ^(b) j f(x)dx, a < £ < b. 


Let us apply this theorem to the Integral 
O+loo 

O^loo 


VTJ 

-^1 


(3.1.24) 

where f(^) Is periodic on the contour 


of Integration, and has a finite meeun value m. 


We write 


f(^) = m + gO). Then 

0<flao 4-00 

O—loo —oo 


(3.1.25) 


It Is now sufficient to ccmslder 
a 


1 g(o+lt) show that It has an upper bound 

Jo Independent of a (we let a -»oo). 


We put 
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o+lt “ o*+i* “ o^+t* • 

We then see that the real part Is monotonle and (3.1.2U) can be 
applied Immediately: 

a 5(a) a 

0 / *** “ 0 / 8 (<»+lt)dt + / g{o+lt)dt. ( 3 . 1 . 26 ) 

Jo Jo 

Since g(o+lt) Is periodic In t with mean value zero, we see that 
the two Integrals on the right are bounded and as a « the second 
term on the right hand side goes to zero so that the upper bound 
for the left hand side is independent of a for a very large. The 
fact that 5 depends on a does not affect the argument because In 
the first Integral of the right hand side, ve can simply choose 
5(a) to make the integral a maximum, to obtain an upper bound. It 
can be shown also that 5 (a) tends to a constant as a 
As for the Imaginary part of 5 ;;^ , ve write 

dt + J dt ( 3 . 1 . 27 ) 

where x Is chosen such that is monotonic in the Interval 

(x, a). The first Integral on the rl^t hand side is obviously 
finite and Is Independent of a, and the second integral may be 
treated as before. 

We have thus shown that If f(3) is bounded and periodic 
on the contour of integration, the integral 


O+loo 


sir / 


tw. 

3 




is bounded. 


0—ieo 


We apply these observations to the integral ( 3 .I. 23 ). There 
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Is a product of two periodic functions in the Integrand* The 
period of one depends on Z* period of the other on H. 

The periods of both are Independent of t). The product will 
be periodic If C is adjusted such that the ratio of the two 
periods is rational* It is therefore possible to state that: 

For every value of $ such that C/H Is a rational n\unber, 
the derivative of Q(l^ Z) with respect to t| la bounded; l*e* for 
such values of Q(]c) has no discontinuities as a function of 

‘x “y’ 

We note that we can approach any point Zt with a g obeying 
the condition expressed In the last statement as close as we 
please* 

3*2 Difficulties of direct use of double Fourier transform of 
Bloch matrix 

Inspection of equation (3*1*10) makes it clear that the 
second Fourier transform may be carried out by use of (3*1*9) 
and one might therefore inquire why we have not simply carried 
out the second transform Instead of working with the quantity 

P)* The reason will now become clear* for we will examine 
the consequence of such a procedure* 

We note that Sondheimer and Wilson's expression Is of the 

form 

c(£i Sa P) » f(H.ri»a) &(£i-£a fi) (g = H 2)* (3.2*1) 
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f l8 independent of p, so that on performing the inveree transform 
we obtain 

p(rira K,) = f(H*Ei * ra) g(ri - r* O. (3.2.2) 

Now the momentum distribution is the diagonal element of 


P(ki 5a) “ I p(Ei Ea) 


a£a 


We put 


g(£i - Ea) =1 Q(k) 


i5*(Ei-£a) 


dk. 


(3.2.3) 


(3.2.4) 


Then 


P(5i 5a) = J <iEi <J£a f(H Ei » Ea) j <55* 
“ J ^(5i + 5 5a + 5) ^(5) 


-i(ki-lc).r, l(^a-k).Ea 


where 


P(ki ka) =.| f e 


-i5ii:i i5a«Ea 

« d£i dra 


«i(5) 

(3.2.5) 

( 3 . 2 . 6 ) 


It can easily be shown that the diagonal element of the double 

Fourier transform of a function of the form P(S.Ei ra) is 

independent of k and k : 

X y 

F(k ^) s C ^(k^) ^(k^). where C is some constant. 

Wo therefore obtain 

P(k ^) a 6(k)C I Q(k^ ky k^) dk^ dk^,; C some constant. 
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[The d-functlon appears since we have taken a double Fourier 
transform and. so obtained the momentum distribution for all 
electrons, and there are an Infinite number of them.] 

The resulting momentum distribution is constant in the 

k and k directions. This is obwiously physically incorrect. 

X y 

We must expect that when magnitude of the magnetic field 
is very small, the true momentum distribution approximates 
closely to that of the electrons in the absence of the magnetic 
field. How are we to explain this peculiarity of the Sondhelmer 
and Wilson solution? 

We think the answer is as follows. Sondhelmer and Wilson 
avoid the troubles earlier workers encountered with boundary 
conditions by working, not with a finite system, but with an 
infinite one. They obtain the solution of the Bloch equation 

0(E' E e) + If (£' E p) “ 0 (3-2.7) 

with boxindary condition 

C(r' r 0) = 6(r^ - r) (3.2.8) 

(H is the Hamiltonian) 

and there can be no doubt that their solution is correct for the 
infinite system. However, supposing Cj|(r' r 8) is the Bloch 
matrix for N-partlcles, the momentum distribution in the limit 
as N-» OB is the inverse Laplace transform of 
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0(JB 3) 1 .. 


T 


N 


I Vs' E P) 


-ik*£* lk.r 

M S* W 


dr dr*. 


(3.2.9) 


Volume of 

■Tstem 


In taking a double Fourier transform of the S-W solution we 
hawe reweraed order of integration and limiting proeess, and 
this appears to he the reason why we obtain an unphysieal answer. 
(It should be noticed that the limiting process here is of a 
different nature to the one we make on equation (3.1.22). Here 
limits of integration wary aa we vary the number of particles). 

In working with the mixed matrix C(k r ^) we have avoided 
these difficulties. Analogous to (3.2.7). the mixed matrix 
obeys the equation 


^ 0<k|E @) ♦ § (SIE .) = 0 

with boundary condition 

ik.r 

C(SlE 0) = e ", 


( 3 . 2 . 10 ) 


(3.2.11) 


ajid it can be readily verified that the quantity defined by 
equation (3.1.10) satisfies equations (3.2.10) and (3.2.11). 


3.3 Relation to Van Hove’s conjecture on use of perturbation 

In concluding this discussion, we shall refer to a conjecture 
by Van Hove (1939) which has only recently come to our attention. 
Van Hove has suggested that whereas it is known that with 
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attractive potentials the zeroth order Fermi distribution cannot 
be used to generate tne correct solution when the interactions 
are switched on, the same might also be true for repulsive 
interactions. 

While we are not able to draw a decisive conclusion about 
the correlation problem from our work above, we can say with 
certainty that in the k and k directiona the momentum 
distribution does not have a discontinuity. However, while we 
note that. If generated by perturbation theory, we can expect 
such a discontinuity to be preserved. It should be pointed out 
that the perturbative treatment leads to a pathological result 
in this case. We do feel however that our ability to carry 
through the calculation of the momentum distribution exactly 
in this case is of considerable theoretical interest. We are 
also of the opinion that the result may have practical significance 
and we are at present planning the numerical calculation of the 
momentum distribution in the presence of a magnetic field. 

Compton effect studies or position auinlhllation experiments in 
the presence of a magnetic field may prove illuminating in the 
future. 
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k. VTlatUnal eeeond-ordT dentlty aatfloes Dnllt trcm 
twe-body ortltals 

We hare eo far dlaeaaaed ehlefly the extreae llalta of 
weak and atrong ooupllng. Unfortunately, the aethoda need 
there are not adaqnate for the range of IntOraedlate coupling, 
and we therefore turn our attention at thla point to acre general 
eonalderatlon of the fora of the aeeond-order density matrix. 

All erldenoe aeeaa to us to point to the fact that thla quantity 
afforda the beat arallable tool for dealing with aany-body 
probleaa In which only two-body forcea operate. 
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4.1 DlseuBBion of nece8«ary and eufflelent oondltlone for 
variational ▼alidlty of Bccond-order denBlty matrlcBB 
In this section we shall fornulate necessary and sufficient 
conditions such that a nany-hody wave function e:q?anded In terms 
of orthonormal two-hody orhltals (rifa)) shall he antisymmetric* 
These In particular provide necessary conditions that the type 
of second order density matrix suggested hy Young and March (1960)* 
viz: 

N^-1) 

r(Si'Ea'| rija) - ^ (E«'Ea') ^^(r,, r,) (4.1.1) 

dal 

(where 2N Is the number of particles) shall he derivable from 
an antisymmetric many body wave fianctlon. In §4.3 it will he 
shown that with the particular form of orhltals chosen hy Young 
and March In method (B) of their paper (which method Is at first 
sight the more promising of the two they propose, since It 
reduces to the correct Hartree-Fock solution in the hl£^ density 
limit) the conditions caimot he satisfied, at least In one 
dimension. 

We will also discuss the question of the sufficiency of the 
conditions on the orthonormal orhltals [^^(Si£a)} Yoz* density 
matrices of the type displayed In equation (4.1.1) to he derivable 
from an antisymmetric wave function. The sufficiency follows 
from an assertlcn hy Bopp (1959). of which Coleman (unpublished 
work; I96I) has questioned the validity and has claimed to have 
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produced a counter exanple* Both the laouna in Bopp'e argument 
and Coleman's eounter-ezaaple will he examined In 

We take a oenplete set of two-hody orthornormal orhitals 
(to he referred to as gemlnals) and e:qpand the 

antisymmetric wave function for 2N particles in terms of them: 


( 4 . 1 . 2 ) 


it •.d. 


also 


^ must he invariant under permutation of and 

N 



The second order density matrix is 


r(Ei'£a'|r,ra) = ^ a^,^ ♦e»*(SiEa)+^(£i£a), 




(4.1.3) 


= H(2R - 1), 

I 

where 


(4.1.4) 


•«.« - scan - 1) ^ 

^a • 

One can see that a^«^ a a^«* - the matrix (a^*^) is Hermltlan; 
it is thus dlagonalisahle and there exis gemlnals [^^(£i£2)i* 
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which we shall call natural geainals, such that 


r(Ei'£a'|r,ra) = ^ ^^(EiSa) (4.1.5) 

I 

Now the antisymmetry requirement 

▼(EiEa«*E2N^ = “ ’(£aS<‘'£2N^ 

Implies that 

y °^ 1 ^a..^J♦^^S 1 Ea) + ♦e^(EaEi)J = 0 (all ,«3. • .-^ 1 ,) 

( 4 . 1 . 6 ) 

and the requirement 

»(EiEaE3E-“*‘) = - ’(EiEaEaE^***) 

implies that 


(all ^ 3 »..»^j|) (4.1.7) 

Equations (4.1.6) and (4.1.7) are necessary and sufficient 
conditions that the expansion in equation (4.1.2) is an allowable 
many-hody ware function. If the expansion is in terms of the 
natural gemlnals, equation (4.1.6) is replaced by 

^^(EiEa) “ " ^^(EaEi) (®11 ^)* (4.1.8) 

Let us now consider a density matrix formed from a finite 
number of gemlnals, 

r(Ei'Ea'lEiEa) - ^ ♦«“(Ei'E8')'t'^(EiEa)* (4.1.9) 
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The eaae orbitals, and no other* or^ogonal to the* appear 
in the wave function, since we hare 


Yj 


Hence 





(Ei£a)- 


.¥ 




( 4 . 1 . 10 ) 


Conditions (4.1.6) and (4.1.7)* where now there is an upper 
limit K on the summations, fora necessary conditions that we can 
choose a matrix of the form (4.1.9). For particular choice of 
a^,^, however, the conditions are not proved to be sufficient. 

In particular. If we take N(2N - 1) emtisymmetrie orbitals and 
prove that an antisymmetric wave function for 2N particles may 
bo expanded wholly in terms of them, we have still not proved 
that these orbitals are natural geminals and so the density 
matrix is of the Young-March form (equation (4.1.1)). 

However, Bopp (1959) states that there is an upper bound 
to the a^: 

0 < 1 (4.1.11) 

If equation (4.1.11) is true, then it is easily seen that the 
second order density matrix for an antisymmetric wave function 
formed from N(2N - 1) geminals must be of the Young-March form. 
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However, equation (4*1*11) is not established beyond doubt, and 
thus will now be discussed. 

4.2 Bopp*s assertion and Coleman*s counter example 

Referring to the eicpanslon of the second order density matrix 
in terms of natural gemlnals, (equation (4.1.5)) it can readily 
be shown that if 

0 <b^<1 , (4.2.1) 

then this also holds for the diagonal elements of the matrix 
in any other representation (that is, equation (4.1.11) is also 
true). 

Let us write the many body wave function as 

* ■ Z Z 

n K 

where Is the set of natural germinals, 

“ 2H 

=TT (l) (K labelling the configuration 

1-3 ^ 

and the are orthonermal and defined by the integral equation 

I (13) #(23) U^^“)(2) d(23) - (1). 

We now antlsymaeiize and so form a normalised determinant 

(3...2H) by which we replaee In equation (4.2.2). 

We now have a prodmet 

#n (3..21f). We antisymmetize this 


( 


I I 
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and again nornallae, obtaining a function 

(1,2.2H). 

Now Bopp glTea an upper Bound for the h^: 

bn < ^ l< ♦ >1* (4.2.3) 

K 

which he replaces by 

b„ < 1. 
n 

It is easily seen that this last equation lamediately follows 
if 

It can be shown that 


♦„( 12 ) - ( 1 ) ( 2 ) 


Im 


where 
and if 


= 0 


all p ^ m 
all q ^ 


“<p - 0 


q^ 


(4.2.5) 


( 4 . 2 . 6 ) 


Thus apart from a normalisation factor Xj^ is equal to 


E V *** (4.2.7) 

V ii ^ 

Consider a second function X^. This is, 

I V E *«■ “<c. — ”*2 j(2»)-('‘.2.8) 

p da 


I 

( I 
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i^re 2^ JJpP 
P 

Now If the products 


iB the entlBynetrlBAtlon operator. 


n “41 \ 


(n) 


ki-3 


#Cj»3 


differ hy more than two functions It Is easy to see that 


< V - 0 


(4.2.9) 


Also, If the two products only differ hy one function then by 
condition (4.2.6) equation (4.2.9) again holds. 

Howeyer, If the products differ by just two functions, say 
^ 0 ^ 3 ( 3 ) 

U^(4) * u^(4) 

then < > can be seen to be proportional to 


® KtiU “ kak* * “ kalu ®/c»K* 


This expression will usually be zero because of condition 
(4.2.6), but not always. Thus equation (4.2.1) does not 
Immediately follow from equation (4.2.3), and In fact Coleman 
claims a counter example to equation (4.2.1) which we now examine. 
First we derlys the expression for the upper bound for b^ which 
Coleman uses, examining the conditions carefully. 

Let be a set of gemlnals, and expand the wave function 
as 

1(12 ... M) - ^ (12) Xn(34 .. M) (M . 2N) (4.2.10) 

n 
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with < ^ > » *„»,» 

aM If 

< X, Xn > - «,„. 

The density matrix is 




the are the natural geminals. 


r(1'2' 112) . ^ 0/ 0^ ♦,( 1 ' 2 ') ♦„(12) <)(,■■)(„> 

an 

Thus if and only if the are natural geminals. 


a 


nn 




1 ). 


We nay then derive an upper bound for a^ = b^^: 

\ = l^nT = ^n^12)Xn(3..M)| f(l2..M) > |“ 

* l< > I* 

(A being the normalised antisymmetrization 
operator for U particles). This is also 

5tn M* < < ^n^ 

The equality will only hold provided 
▼ a XqI . 

Now Coleman takes a 4-particle wave function of the form 

▼(1234) B (iA|^( 12)^(34)1 . ((i some normalisation constant) 

and shows that with a suitable choice of orbital ^(12), 

< ^(12) ^(34 )|a|^(12 ) ^(34) > 

can be made arbitrarily close to ^/3, where Bopp's upper bound is ^/6, 
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However, we have remarked that the hound for a^ only holds if 
la a natxiral geminal, and this la what Coleman must assume 
for his counter example to he valid. However, he fails to prove 
that the orhitals he chooses are natural gemlnals, and without 
such a proof his counter example cannot he aecepted. 

We therefore conclude hy remarking that neither condition 
(U.2.1) nor Coleman's counter example is hevond douht, so that we 
cannot provide conditions which we definitely know to he sufficient 
for the variational validity of suitrlees of the form suggested hy 
Young and March. In view of the suggestive and promising nature 
of the form expressed in equation (U.1.1), further work in this 
direction would obviously he valuable. 


4.3 Two Body Bloch Orhitals 

To introduce the Wlgner lattice of the low density limit we 
will take the two body functions to he 


8u»Sn 


(4.3.1) 


where the are reciprocal lattice vectors. 

When the nvunher of particles becomes very large, the many body 
wave function T(£i£a..£n) approximated to hy 


j ®(j5i5a|l5sS4| .. 

% 


ja (£i£a) .. 


(4.3.2) 
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1 


(n^ l8 the ToluBe of the Brlllouin sone). 

How we Bust entlsyBBetrlze in £i end £•, etc. We can do this 
hy mahlng 

C(]si)cal...) - - C(]c«ki I...) (U.3.3) 

(here and helow the unepeelfled Indices are assuBed to take any 
value, provided we inpose the oondltioa 

▼(j5i|ka) » v(ka|ki) (4.3.U) 

Equation (4*1*7) hecomes 

Yj j C()Si)Salfe3fe*l *Ov(lSi+KBll5a+l5n^^(S»+V^*5v) * 

e “e “e **0 

-l(Xi+Ka),ri -l(jCa+5b).ra -I(ft+Ka).r3 
Multiplying hoth sides hy e e e 

'S'* 

a e ^ and integrating over all space, we obtain 

6-functlons, giving 

C(SiSal £3«41.. )v(£, +5,^1 Sa+5b)^(«»+5al £*+8p) 

- - C(£i£3|(Cafi4| ••)T(Si+6jg3+K^)v(5a+5blS*+^) (4.3.5) 

Now for one-dlaenslon the functions proposed hy Young and March 
obey equations (4.3.1) and (4.3.4) if we force periodicity in ki 
and ka (although this is not essential to the arguaent, only 
affecting the precise natxire of the labelling of the v(ki ka) 


I 

I I 
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with and Kjj), for, putting y(x) ■ p(x) + x, p a periodic 
function, the ahoye author's functions (equation (5*2) of their 
paper) nay he rewritten 


with 


♦k, 


(xi-xa 



l(ki-ka)P(Xi-Xa) ikiXi 
e e 


ikaXa 

e 


(4.3.6) 


It will he seen that this is of such a form that 

v(ki+Kj^|ka+Kj,) » 0, ^ 0. 

Hence from equation (4.3.5) 

C(kikalkak4|..)y(k,+KJka-K^)y(k3+K^|k4-K^) - 0 

unless K. « K„. 

a a 

Now if we are to obtain a second order density matrix of form 
(4.1.1), then we must he able to find for a given k, ka(ki#ka) 
a value of and values of ka ka ... k2jj such that 

C(kika|kak*|..)v(ka+Kjk*-Kjjj) # 0, 


But than 

v(ki+K^|ka-K^) ■ 0 unless 

This shows to he singly a product of two plane waves, 

which result must hold for all ki, ka. We have thus shown that 
method (B) of the Young>4lareh scheme in one dimension is only 
varlatlonally valid when it reduces to the Eartree-Foek result 
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(th« wave funetion than balng almplj a dateralnant of plana 
aayaa). 

It muat ha ranarkad that Yeung and March wara unabla to prora 
the satisfaction of the Pauli conditions on the first order 
density matrix in thalr method (B)» HoaaTar, it is known that 
satisfaction of the Pauli conditions does not snsuirs a 
varlatlonally valid scheme, and it would appear that we muat 
ensure the ezlstance of a wave function directly. Nevertheless, 
the above example shows that it may not be necessary to explicitly 
build up the wave funetion to ezaualne the validity of particular 
trial forme of the two«-body orbitals of equation (4.1.1). 


44 



5* Localised orbital method (with M. Darby) 

We turn now to discuss an explicit nethod of setting up a 
trial seeond-order density matrix. No really accurate form 
for the whole range of electron densities has, as yet, been 
found, but the method outlined here should be valid for densities 
irtxleh are not too high. 


5»1 first and seeoBd-srder matrices 

The energy of the uniform system of any density may be 
written as 


s 


- jir / [y T(E'£)] ^ / 4E' I SE 


t2P(g*e) - 1] 
lE'- El 


Where the energy Is in atomic units, and 


X 



9 


(5.1.1) 


r^ being the mean Interparticle spacing, and r(£^£) snd the 
pair function P(£*£) are for a gas of uniform density. We 
take for the paramagnetic gas 


2r(£i£a' £i£a) - ^ ' «*(£a'- £i-Sn? (5.1.2) 
8n 

The 2^ are vectors of a body-centred lattice, and the prime on 
the summation sign Indicates that ■ 0 Is exiluded. a(£) Is 
a Wannler function for a simple cubic lattice (the b.c. lattice 
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'being aadc up of two of the eluple ou'bie lattloee) whloh hue 
nodes at all the lattice polnta of the alaple ouiblo lattice* 

We fopa the first order density matrix froa 

Y(£a'£a) J ^(li £a' ra)dri. (5*1*3) 

We write 

s(e) “ ~ (5.1.4) 

Is the Brlllouln zone for the simple cuhlc lattice (the rectors 
of which we shall denote hy S^)* Inserting (5.1.4) In (5.1*3) 
and letting N -» « we obtain 

/ X r . / X * ^)5*^£i“£a) 

y(£i £a) = J l▼(J5)|“ e dk. 

The kinetic energy is therefore given by 

“ ^ / ^*Y(r'r)j.,^ = j ar j |y(k) |» k» dk, 

so that kinetic energy per particle is 

“ -ji*- ^ at 

which Baybataarrltten 

T » - 2 p- j a* (£) V* a(r) dr (5*1*5) 

It may be noted that we antoaatlcally satisfy the Pauli condition 
that the occupation numbers are less than one* 

0 < a(k) < 1, 


I 

I I 
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f«r we can fern a Blooh fonctlen 

>)5-Sn 


“ Yj *^S-8a)® * ^ ▼(JE+Sn) 

^ Sn 


i(5+5n)* 2 


and this nonialiaea to unity over the unit cell: 

1 - / «£■“ ^l''(!S+Sn)l** 

Kn 


Henoe 


T(j5) < ererywhere. 

Tor the potential energy, since P(£'£) » P)r'>r), we obtain 


per particle 


V ■ 


i / 


[2P(r) - 13 




>ow 


2P(£) • X! * '•“ 

Bn 

We therefore hare 

I^Ce)!* 


[ 17 

Bn 




Ir<„i 


( 5 . 1 . 6 ) 


Hence the potential energy la Juat that of the orbital a(£). 

centred at the origin In the fields of point charges, one at 

f ^ 

ewery other lattice alte (the tern j __ will be Ignored below. 
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cancelled the aelf energy of the unifore hackgroond change). 
If we fora the periodic potential 

Bn 

we aay write the total energy per partlele aa 

/ / ^ ds. 

(5.1.8) 

where we have taken the mean value of V(£) aa zero. Thla la 
of no consequence when we seek Euler equations with particular 
trial forma of a(r) such as the one we now examine. 


5.2 The Macke Tranafcreed Wanni.it. 

In order to generalize the scheae outlined above. we 
replace the functions a(r) by Macke-traneforaed eii^le cubic 
Wannier functions, given by:- 

(XiYiZi)* a^(X) ay(Y) a^(Z) - a(R) (5.2.1) 

Here X = X(x), Y = Y(y), Z = Z(z), emd we restrict these 
tranaforaations so that Xi(x+Sj^) = X,(x); Xi(x) » X,(-x). 
Further Xi » dX/dx; Xa = d*X/dx*, etc., and the simple cubic 
Wannier function is defined by 


S(x) 




dk 


11 and flg are respectively the voluaes of the unit cell for the 
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elBple cubic lattice, which has yectors and the Brillouin 
zone for the same lattice, and are thus related through:- 

flflj ■ 8**. 

It la conTenient to quote here certain sunaation rules 
for the siaple cubic Wannler functions. Defining the length 
b by Dq ■ 8b*, we hawe 


^ “*(E^) «<E-Sb) - Vn. 

Sn 

^ to - Yj ^ “(E“§n^ " ° 

2n 2n 

Y ^ a(E-2xi) * " ^ to «“(E-2n^ to 

Sn §n 


These will be needed in obtaining Euler equations, but before 
we do this, we atust see what energy expression (3*1 *8) becoaes 
when we replace a(£) by the transforaed function (3>2.1). 

The kinetic energy tera consists of an integral ower all 
space, the x-eoaponent of the Integrand of which is just, 

where • da^to, etc. Siallar teras hold for the y and z 
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components* Now since X(x) Is periodic In the lattice, we nay 
replace the Integral over all space hy a sum of Integrals oyer 
one zone, and then employing the eummatlon rules (5.2.2) we 
Obtain for the kinetic energy per partlole:- 

jibr / [t - Sr - S 


(5.2.3) 


In order to obtain an Euler equation we require the energy to 
be In the form of an integral over all space, and thus we 
Introduce the first of the summation rules (3*2.2) Into (5*2.3)* 
obtaining:*- 




o*(£)o(r)dr + similar terms for 
~ ~ y and z components. 


All space 

Replacing a(£) by the transformed functions (5.2.1) In the 
Potential energy part of (5*1.8), we obtain 


/ ar [V(r) - 1 ] Y,Z,X, a»(g) a(g) 
All space 


* O ^ 

All space 

Here V'(£) Is a localised potential. Again using the first 
summation rule of (5*2*2) on the first term we obtain finally 
for the potential energy per particle 
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(5.2.5) 


^ J V(r)XiY,Zi a*(£)a(r)ap - ^ J J X,Y,Z, a*(g)a(R)dr 
All space All apace 

The total energy per particle, using Maeke transformed 
simple cuhlc Wannler functions In this scheme. Is thus given 
by the sum of (5.2.U) and (5.2.5). We proceed now to obtain 
the Buler equations for the transformation functions X,Y,Z. 
Functions satisfying these equations will give the lowest bound 
of the energy as given by (5.2.U) and (5.2.5). 


5.3 The Baler Bcuatlons 

If we write the energy of our system In the form 

Cm I I(X,Y,Z,X,,Yi,Xi,Xa,Ya,Za)dxdyds 
All space 

then the fonctlens X,Y,Z making 6 a nlnlmun satisfy the Buler 
equations for oursdtas. The variations In X,Y,Z must be 
considered separately, and In partleular, for the variations In 
X we have:- 

00 OO 00 

■□■JdxJdyJdsdj^ I 

•00 * 00 —oo 

Assuming I to be sero at the limits x > - «, we may write 
00 00 

J 4y J ds * 3^1 » 0 (5.3.1) 
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Ab X -» X -f 6X; Xt Xi + 6Xt j Xa Xa + 6Xat ■thus: — 

■ S “ * # * Ifc **• 

Further 

fit ^ ^^ [“jT -fhit “^ 

By our previous aesvunptlons the first tern on the R.H.S. is zero . 
Similarly 

f or = I ^ ^ 6X d£ 

Nov since SX is arhltary, the Integrand with respect to X must 
vanish for all X, (hut only X) so that (5.3.1) becomes 

l^f^[§--A ^ * 8? ^] = ° ‘5.3.2) 

similar expressions hold for variations with respect to Y and Z, hut 
here we will confine our attention to those with respect to X, 
the Euler equation for Y emd Z being obtained by rthawgi ng the 
variables in a cyclic way. 

In evaluating (5.3.2) it is convenient to consider the 
kinetic and potential energy terms separately. For the kinetic 
energy part, I is the integrand of (5.2.4). Bvaltxating the 
derivatives and combining terms ve have for the kinetic enez*gy 
contribution to the Exiler equation:- 
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00 00 

^ J ay J dz YiZ, [a«(£)[^ [- 5^7 


Xa* 

W7 


Xa* 

- 5x7^ 



«• 00 —00 

3^ X ^ 

t 2o'(i)a(£) [- 5S7 + S7T 




Xa“ 1 ,> 

♦ ST' *X 



♦ ^ ( 2 “'(i:)o.(e))[-^]] 

The ixiflnite integrals may now he written as sums of 

integrals over the sides of the simple cuhic unit cell. Since 

we have only two integrations we obtain a sun where 

8 .8 are the y and z components of 8 respectively. In order 
y z •*** 

to use the summation rules (5.2.2) it is convenient to Introduce 
a further summation, over 8 ^, hut it should he noted that this 
now has to he carried throughout the working. Using the 
summation rules we obtain 

"v]]/ ■»'«»/ (5.3.3) 

-a -a 


where a is defined by H s 8 a^. 

Using expression (5.2.5) we find that the potential energy 
term contributes 


OD 00 

it I*'I If* [»*(») ^ (i) -«*(£) & T(B' 


- T(b) 2«(*) «'(S)] 


For the terms involving V(£), the infinite Integrale can again 
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tf replaced tj sue of integraie orer the sides of the unit 

cell, and to allow use of (5*2«2), and for consistency with 

(3*3«3) ve SUM ower Thus we ohtalnt- 

a a M 00 

- sir / ^ / ^ If-ZiVCE) Y,f ^ Y,Z,a>(Ii-3jP^ 

-4i -a - — 00 - 00 

°x 

(5.3.4) 


The Euler eg]aation for X is now obtained by putting the sum Of 
(5.3.3) end (5*3.4) equal to zero. Before doing this we note 
that if we choose the density of our system, Y 


T 



h»k « 

5 ? 




equal to unity, then we flnd:- 
n a 2.; a 


= 2"*“; b * *2"^ 


Then we have;- 
a a 

Yidz ' 

-a 


j Yidz I Zi dx = 4a* = 2*’ 

-a -a 


Finally, since XiYiZisre assumed to hare the same functional 
fora, F say, we may generalize the above terms for Xi and obtain 
for our transformation function F the Euler equation:- 

2 ^ d r *'»(*) Pa*(x) 1 -| 1 f * 

s L“ y 7 (x) + siTTx) + 2 *itPi»(x) J “ ijc J ^ Pi(y)Pi(z)v(£) 

y Pi (y)Pi aMP(x) - s } - 0 

i [(x-S^)*+y*+z*]“ * 

(5.3.5) 


Q *00 • 00 

®x 
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This equation glYes the heat transformation functions F* 

It nay he ohserved that in the limits of high and low densities; 
F(z) - z gives the plane wave result in the limit X -» 0, and 
Fi(x) ■ ^ 6(z - is a solution in the limit X -» «• 

1 Work is proceeding on the detailed properties 

of the solution of these Euler equations* 


55 



6. Non-\mlform electron gaa (with 8. Baapanthar ) 

So far, we hare eonaidered the Influence of eozrelatlon 
effeete In a uniform electron gaa. Unfortunately, thla 
Soamerfeld model neglecte certain featuree which are of central 
Importance In the theory of metala, and In particular areragee 
out the periodic potential due to the lone to a constant walue. 

Very recently, Bellemana and De Leaner ( 196 I) hawe reported 
a many-hody approach to the non-uniform gas problem. Their 
results, which were simply quoted without proof, have apparently 
been obtained by expanding the grand partition function of the 
system in powers of the coupling parameter between electrons *t><< 
the posltlTe (point) ions and between pairs of electrons. By 
summing infinite series of divergent terms, they were able to 
obtain finite restilts for the energy per particle. We show 
below that, by finding Dirac's density matrix in a self- 
consistent framework, we obtain a resiat which is closely 
connected to their final energy formula. Indeed, this may 
be obtained almost at once for our theory, if we add the 
Oellmann-Brueckner correlation energy for a high density uniform 
gas to the energy calculated from our approach. 
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6«1 Bneggy and Dlrae dwl^y matrix of a non-uniform electron gaa 
We ahall eonslder the problem of electrons moving in a 
periodic lattice of protena, in order to eircumvent at present 
the difficulties of introducing Hartree fields associated with 
innar eleetrmna, and for the eoaduotion electrons we shall 
adept essamtially a Hartree description. We first write down 
the diagonal element pC?) of Dirac'a density matrix, **^^0 firat 
order in the Hartree periodic potential V(^), uaing the density 
matrix perturhatien theory of March and Murray (1^61). If 
ia the magnitude of the wave-vector at the Fermi surf see, and if 
we use atomie units throughout, we have then 

I- / d?. —: 

I? 

Ji(p) • (dim p - p coe p)/p** We see from (6.1«<1) 
that the eonstant density pe s lc^V3<* Im modulated hy the 
periodic density eerreetion, which we denote hy pi • We now 
introduce the Fourier components of the Hartree potential 

Y(?) throng 


P(?) 


V 2k • . Y(?,) da2k^l? -^il) 

- ^ - 


( 6 . 1 . 1 ) 


▼(?) - ^ 0 ^^^“*^^ 


( 6 . 1 . 2 ) 


where the are reciprocal-lattice vectors, and we impose self- 
consistency through the Poisson equation, 

»*y(?) - iw 7 •(? -1^) - iwpff) 

«n 
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- 4*p(?), 


(6.1.3) 


¥ 


l3! J? 
& 


wlxor* th« danote rector a in the direct lattice and Q. la the 
n 

volvuae of the unit cell. Hence, auhatitutlng (6.1.1) and (6.1.2) 
In (6.1.3) and ualng the reault that 


2k^* f 

I 




l#Cj^.? di(2kfr) 


* V «n 1-*3 


we find 




0 , 




0 . 


( 6 . 1 . 4 ) 


(6.1,5) 


( 6 . 1 . 6 ) 


Thia reault, comhlned with (6.1.1) and (6.1.2), deflnea the 
self-oonalstent electron density In this approxiaatlon. 

We next write down the potential energy U in the Hartree 
approximation. If and denote the potentials due to the 
nuclei and the electronic charge diatrihution reapectlrely, than 
V B Vjj -h and it may readily he shown from the full Hamiltonian 
that the potential energy is given hy 

^ * I ^ P ^ I ^ P (6.1.7) 

where is the Coulomh Interaction energy of the unscreened 
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protons; the second term Is the electron-nuclear interaction 
energy and the third the electron-electron potential energy. 
Separating off the Madelung energy of the proton lattice in a 
uniform baekgreund of electrons with density po> ve may rewrite 
(6.1«7) in the form 

D = U^adelung + i ^ Pi <?) V(?) + i J d? p, (?) Vjj(?), (6.1.8) 

where we have used the fact that j d? pi (?) =0. To obtain the 
total energy B we must add the kinetic energy contribution, which 
in the density-matrix perturbation-theory has been given by 
Corless and March (1961) as 

T ■ jj ^1 V(?i ) V(?a) 


Ji (2ky ?a ■ ) 

* (2kj. ?a - ?, )* ^ 


o(v»), 


(6.1.9) 


where cue first term is the usual Fermi energy for free electrons. 
Rewriting (6.1.9) in terms of pi defined from (6.1.1), it is 
immediately seen that the second term in (6.1.9) cancels the 
second tern in (6.1.8) when we form B s T TJ. Thus we find 


* * ’^Ferml * ^Madelung + i J ^ Pi (^) (6.1.10) 

The last term ean be evaluated by using the Fourier-aerles form 
of Pi and the fact that V-(?) is a sum of Coulomb potentials over 
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the lattice sites. The contrlhutlon from this term to the 
energy per particle may then he readily obtained as 


V(2 ) 


(6.U11) 


Adding the usual exchange and correlation energies to the 
Hartree energy derived above, «e obtain the resxilt of Bellemans 
and De Leaner if the Fourier comp^ents V(/c^) In (6.1.11) are 
replaced by the unscreened Coulomb values 4A^’/3 ic . Our 

result appears to have more direct physical significance in that 
screening is Included from the outset. 
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6.2 Discussion of natural apln orbitals 


A complete set of Bloch functions will he defined 

hy the expressions (for an infinite crystal) 

“ Z(6.2.1) 

Sn 

/ “(S,. 8n>* <®****> 

All space 

It follows that 


Now 


/ I V<e>l * ^ ' 


the normalisation is Independent 
of Jj. 


/ ♦*/ ‘(E) (E)^ = ^ / .l‘(E^) .j(E- 

All space Bm»Sn 


■H^Je 

-n' 


- ^(j5s“Si).5n 

* I, ' 

8n 

/ ♦l5,‘^(®)V^'*)*® " ”b *11 Xi *'*’ ■>•**»)• 

All space 

t(e»E*)» defined for a finite crystal, will obey 
periodic boundary conditions analogous to those laposed on the 
one particle Bloch functions. These fora a complete set so that 
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y(£' ,£) Miy be expudad in t«rs« of thorn over the domain of the 
oryetal. We ehall treat the infinite eryatal eaae* We write 


y(e',e) - y / <1)51 / eij()5i.)5t)*v,“*(r')V ^(C). 

li \ \ 


(6.2.4) 


To obtain the natural orbitals we auat diagonalize C.^Cki.ka). 

We shall show it to be already diagonal in j^. 

Equation (6.2.4) only requlrea C^j(ki.ka) to be defined 
when and jj^a within the Brlllouin zone. However, we ahall 
for mathematical convenience take it doubly periodic in the sense 


Cij(l5i+l..)5a+5ti^ “ Cij(lj,,ka;. (6.2.5) 

Our results will follow from the imposition of the 
periodicity condition expressing the arbitrariness of origin in 
the crystal; 

Y(E'-«n»2“8n^ “ ^^S'»S)‘ (6.2.6) 

We define the mixed matrix 

t(s.e) - j y(e'.E) d£'. (6.2.7) 

All space 

Thus 


Y(eS.S+Sn) - / Y(£',E+Sn)® ‘^*^ ^E' = / y(E'+Sn»E*Sn)®^"^“ 




y(s»E)«^~‘^ » Y(«,E+8n)* 


( 6 . 2 . 8 ) 
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'A 


Now sine* Y(<fV) is ahown to be an eigenfunction, elgenTalue 
• of the traaslation operator 1^, 


Y(£f£) - ^ Cj(£) 
1 


(6.2.9) 


Let us BOW write, 


Sn 




IB eqjaatlon (6.2.3). We see that 

Y(£.r) - 8** as Cij(£.jc) V^*(<) 

id % 

er 

^Oj(s) Vcs) - as 0jj(5,S)Vj»(s)»ji(j). 

J Ij % 

Nultiply ttureus^ by ♦|j5^*(£) “8 integrate over all space. 

(6.2.3) shove that 

^ »(S«-S*Sn) - 8»* ^ J aj Cij(5.S) y '>(>'o-S*Sb)Vj“(s) 

^ ^ *1 ”® Sa 

OF OjCs) ^ »(s-s*li.) • ^ eij<S.S) V<s)- 

I. i 

IntratiBg over the Brllloaln zewe with respect to | and patting 

Si 
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I 


OjCg) - 8*» ^ Ojj(£) Vi"(£), 


( 6 , 2 . 10 ) 


It should he noted that equation (6.2.3) Inpllea the 
periodicity of Cj^j(/g): 

°i3<£*5.> - °i3<£>- ' 

Our expression for Y(<t£) nay now he rewritten as 

y(£»E) ■ ^ ^1 


( 6 . 2 . 11 ) 


Hence 


y(s'»E) » / y(S.E)« ^ ^ ^‘YjI 


-lic.r' 

MM. 

e Ok 


u “b 

1 *0 • 

r(£',s)‘Y,l =13(£V"‘s')V<£’'*S- (6.2.12) 

13 ^ 

Let us consider the effect of applying a unitary 
transformstion 


where 


Y, ®id (£).^ 

*lk “(*) •!.(*) ■ »»<. independent of k 


and also define 
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n m 

Than 


% 


(£)^ 


K^ Ij) % 
mn 


■ E / <=13<!s)*3n‘i. 


mn 


tj“(E' )1 's"(e)4S 


'1/ 

id “b 

We now prove that the set l^jg^(g)l conelste of Bloch functions 
hy forming 

and demonstrating that 

/ “<*(E'«,)«,(E-%,)4£ - % ♦(S,,8„) 

All space 
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We haTe 


a ^ n Ry % 

Hence 


%• / “/(£-8,>“p(i;-Bn)^ "EE/ “’ / 


«s.‘)S'‘^'Sp'ds*) 


Ika .R 


/ *n*(2-Bn-Sm^®n^S"Sv“Sn>^ 


v~' r 


LLI 


dk 


Sv 




ki/ 

% 




Or 


/ “«*(E-SB)ae(j;-Sn)‘^ - / dk, e 


lki(R„-R„) 




m 


kl. 


dki e 


i)Si (Sn-Sm^ 






The asaertion is proved* 

Now we note that because of the essential equivalence of 
k and % + Kj^, a%j,(k) can he taken ae periodic in k, 

«e see that is a set of Bloch functions for all 

unitary transformations, including that making Cjj(k) diagonal. 
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Thus 


y(s',e) 

1 ^ 

and the natural orhitala are Bloch functions. 

To Buminarlzey we can write 

Y(s'.E)=yf <tti f <*^8* 0^3(151 »ISa)t)5/*(E')+)5a^(E) 

u\ \ 

where la a complete set of Bloch functions. 

A 

C 44 (ki,ka) Is shown to he diagonal with respect to ki and 
hy consideration of the mixed matrix Y(<c»3f) which Is such that 

Y(5,E+Bn) « y(s,s). 


We show that if we apply unitary transformations (aj^^(^) to the 
(Cij(JS)) we obtain the set » 2 
te be Bloeh fonotlona. 

The occupation numbers of the natural orbitals we obtained 
by applying unitary tranafomatlons dlagonallslng the (C^j(lc)). 
The natural orbitals, obtainable by the same unitary 
tranaformatlon auat therefore be Bloch functions. 
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